Field and intensity correlations in amplifying random media 
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We study local and nonlocal correlations of light transmitted through active random media. 
The conventional approach results in divergence of ensemble averaged correlation functions due 
to existence of lasing realizations. We introduce conditional average for correlation functions by 
omitting the divergent realizations. Our numerical simulation reveals that amplification does not 
affect local spatial correlation. The nonlocal intensity correlations are strongly magnified due to 
selective enhancement of the contributions from long propagation paths. We also show that by 
increasing gain, the average mode linewidth can be made smaller than the average mode spacing. 
This implies that light transport through a diffusive random system with gain could exhibit some 
similarities to that through a localized passive system, owing to dominant influence of the resonant 
modes with narrow width. 



Without a counter-part in the electronic system, the 
coherent amplification of light adds a new dimension to 
the fundamental study of mesoscopic wave transport. An 
inherent quantum/ wave signature of mesoscopic trans- 
port is nonlocal intensity correlation 0- 0- E| ■ wn i cn re ~ 
fleets the closeness to Anderson localization transition . 
Light transport in an amplifying random medium expe- 
riences enhanced contribution from long paths @, that 
should have a profound effect on the nonlocal intensity 
correlation. 

Due to formal similarity, it is tempting to treat a ran- 
dom system with gain as if it had "negative absorption" 
0, and directly adopt the results obtained for an ab- 
sorbing system. Such simplistic approach to correlation 
functions (CFs) is fundamentally flawed. Theoretically, 
the spatial and spectral CFs are obtained by average over 
an infinite number of random realizations. Among them, 
there exist rare configurations containing more localized 
modes that could lase in the presence of gain. Light in- 
tensity in the lasing configurations diverges, so do the 
ensemble-averaged CFs. Experimentally, the divergence 
of laser intensity is prevented by gain saturation. Nev- 
ertheless, the lasing configurations have much higher in- 
tensity than the non-lasing ones, thus they dominate the 
CFs. The width of spectral CFs is simply equal to the 
lasing mode linewidth, while the spatial CFs only reflect 
the spatial extent of the lasing modes. This is in contrast 
to the "negative absorption" model, which does not con- 
tain the divergent contribution of the rare events In 
order to obtain the CFs that reflect light transport in 
amplifying random media, we introduce the conditional 
average over all non-lasing configurations (...) — > (...) c . 
Such replacement, together with the fact that the fraction 
of lasing configurations varies with the amount of gain, 
make any analytical derivation challenging. Numerical 
simulations turned out to be a fruitful alternative. 

In this paper, based on numerical simulations, we 
present a phenomenological analysis of the local and 



nonlocal correlations of light transmitted through ac- 
tive random media. The systems under consideration 
are in the diffusive regime, but not too far from the lo- 
calization threshold. We show that "negative absorp- 
tion" formulae give a good fit to the conditional CFs 
only at low gain, with a decrease of the dimensionless 
conductance g this range of applicability is further re- 
duced. For high gain, even after discarding the con- 
tributions of lasing configurations, the long-range cor- 
relations are significantly stronger than the prediction 
of "negative absorption" theory, especially for the sys- 
tems closer to the localization threshold. This is sur- 
prising in the sense that removal of all lasing configura- 
tions, which are dominated by more localized modes and 
have stronger nonlocal correlations, should have weak- 
ened the long-range correlations averaged over the rest 
non-lasing (less localized) configurations. Moreover, the 
average mode linewidth 8v, found from the width of con- 
ditional spectral CF, does not exhibit any widening com- 
pared to the "negative absorption" expression. This is 
unexpected because exclusion of the narrowest modes 
(which have lased) should have overestimated the "av- 
erage" mode linewidth. Therefore, the enhancement of 
long-range correlations and narrowing of spectral corre- 
lation width caused by coherent amplification exceed the 
expectations from "negative absorption" model. It re- 
veals the absence of duality between gain and absorp- 
tion. We also calculate the effective Thouless number 
8 = 8v/du, where dv is the average mode spacing. In 
the absence of gain, the onset of localization is marked 
by 5 = 1. We show that for diffusive systems, as gain 
increases, 8 decreases monotonically to below 1 before 
reaching the diffusive lasing threshold predicted by the 
"negative absorption" theory that ignores fluctuations. 
This is an intriguing result, which seems to imply that 
transport in a diffusive system with gain could exhibit 
some similarities to that through a localized passive sys- 
tem. 
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In our numerical simulation, we consider 2D random 
medium in a waveguide geometry, shown in the inset of 
Fig. QJi. It consists of a metallic waveguide filled with 
circular dielectric scatterers of refractive index n = 2. 
Our numerical method for calculation of CFs in passive 
systems has been described elsewhere . Physically, our 
system is quasi one-dimensional, and the transition from 
diffusion to localization can be realized by increasing the 
length L of the random medium. To demonstrate the in- 
dependence of our results on microscopic structure of the 
system, we varied both the filling fraction of scatterers 
and length of the random medium to obtain samples with 
g from 2.2 to 9. The effect of absorption or gain (inside 
the scatterers) is treated by classical Lorenzian model0 
with positive or negative conductivity. The advantage of 
our numerical model is the ability to introduce spatially 
uniform gain as well as to separate coherent amplification 
of an input signal from spontaneous emission of the ac- 
tive medium. In the numerical experiment, a short pulse 
was launched via a point source at the input end of the 
waveguide. Fourier transform of the electromagnetic field 
at the output end gave the CFs for field [ic| C E (Ar, Av) 
and intensity C(Ar, Av). In the presence of gain, long 
after the short excitation pulse, the electromagnetic field 
decays with time in the non-lasing realizations, while it 
keeps increasing in the lasing ones. We excluded the las- 
ing realizations from the ensemble average for CFs. 

Based on the pairing of incoming and outgoing chan- 
nels, three ^lj contributions to intensity CF have been 
identified 0: a local (short-range) C\ w | | 2 , and two 
nonlocal C2 (long-range) and C3 (infinite-range) ones. 
For diffusive transport 5 > 1 in a wave-guide geometry 
C\ ~ 1, C2 ~ 1/g and C3 ~ l/g 2 , making the values 
of C2 and C3 small ; 3|. The nonlocal terms are brought 
about by long propagation paths which are most sensi- 
tive to the effect of amplification. Despite of their small 
values in the passive systems, we observe a dramatic en- 
hancement of the nonlocal correlations by coherent am- 
plification. 

Spatial field CF in 3D bulk was originally found by 
Shapiro in Ref. • Later Eliyahu et al\l2\ calculated it 
at the output surface of 3D random medium. Similarly, 
we derived the corresponding expression in 2D case 8] : 



C E (Ar) 



n(z b /l)J (kAr) + 2sin(fcAr)/fcAr 
nz b /l + 2 ' 



(1) 



where k — 2n/X, I is the transport mean free path, Jo 
is Bessel function of zero order, and the extrapolation 
length Zb ~ / accounts for boundary effects. The imag- 
inary part of Ce{At) should vanish due to isotropy[13j, 
which is confirmed by our calculation where its value is 
less than 10 -3 . The real part of Ce{At) are found un- 
changed in the presence of gain or absorption, as shown in 
Fig. nji for a system of g = 2.2. Eq. Q gives excellent fit 
for passive, absorbing, and amplifying systems with the 
same value of z b /l. Physically, this invariance can be ex- 
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FIG. 1: (a) Real part of Ce{At), (b) real (open symbols, 
solid line) and imaginary (solid symbols, dashed line) parts of 
Ce(Av). Circles, squares, and triangles represent numerical 



data for passive (T abs/amp = 00), absorbing {r ab 



,J ' amp ) 5 



and amplifying (r amp = 5r^ p ) systems, respectively, r^mp 
is the critical amplification time constant at the diffusive las- 
ing threshold. L=0.4 m, W = 0.2 m, I = 0.18 cm. Curves 
represent theoretical fit with z b /l = 0.8. The inset is a sketch 
of the numerical experiment. 



plained by the local nature of Ce{At). Spatial field cor- 
relation contains information that comes from the length 
scale of mean free path. I is always shorter than the 
ballistic gain length l g : l g /l > (2n^ f /ir 2 ) ■ (L/l) 2 > 1, 
because the system is below the diffusive lasing threshold 
L < irlamp- In the above expressions, the amplification 
length lamp — \/DT amp , where D is the diffusion coef- 



ficient and Tamp = lg/c, the effective index of refraction 
nfff — c/ve, where ve is the energy transport velocity. 
Since amplification occurs on the scale much longer than 
I, it has negligible effect on short-range transport and 
local spatial correlations. 

The spectral field CF Ce{Av) contains an important 
dynamical parameter of transport - the diffusion coef- 
ficient D = ve 1/2. The spectral correlation width 8v 
is defined as the width at half maximum of |Cb(A^)| 2 
divided by a numerical factor 1.46. In a passive sys- 
tem, 8v is equal to the average mode linewidth D/L' 2 , 
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where L' = L + 2z b . Since ve can be determined sepa- 
rately through calculation of energy distribution between 
air and dielectric scatterers[^, the transport mean free 
path was found by fitting of the real and imaginary parts 
of C_b(Ai/) [Fig. \I]p]. The value of I allowed us to de- 
termine g = (ir/2)n^) f Nl/L', where N = 2W/X is the 
number of propagating modes in the waveguide, and W 
is its width. In the presence of absorption|14j . the numer- 
ically calculated C_e(Ai/) fits well the expression derived 
in Refs. . For the case of amplification, we obtained 
the "negative absorption" formula by making the substi- 
tution -T abs — * T amv : 



Ce(Au) 



sinh(q a) sm(L'/l amp ) 
smh(q L') sm(a/l amp ) ' 



(2) 



where q = 7 + -*7_, 7± = (^l/Cp + ^Tl/'Lp) A 

(3 = -i i /2ttAj/JD, and a ~ I is randomization length. By 
fitting Cb(Ai/) with Eq.©, we obtained 8v, which is 
plotted in Fig. Elfor systems of g — 4.4 and 9.0. The 
narrowing of spectral correlation width by gain is due to 
partial compensation of light leakage through the system 
boundaries. Absorption, on the contrary, introduces an 
additional loss mechanism, that leads to an increase of 
dv. For both amplifying and absorbing media, the calcu- 
lated 5v agrees well with the diffusion prediction. This 
agreement is surprising for the case of amplification. The 
"negative absorption" theory neglects the fluctuation of 
lasing threshold, and assumes the spectral width of all 
modes decreases with gain uniformly. However, the width 
7 of resonant modes has a distribution P(7), schemati- 
cally plotted in the inset of Fig. For a given amount of 
gain, the modes with small 7 in the tail (fii) of P{"f) lase, 
and they are excluded from the ensemble average. Such 
selective elimination of the narrowest modes should have 
led to an overestimation of 5v. The absence of devia- 
tion from Eq. indicates amplification not only reduces 
the width of all non-lasing modes, but also enhances the 
weight of the modes with narrower-than-average width 
(in Q2) in the averaging. 

In Fig. |2 the diminishing correlation width in am- 
plifying system signifies approach to lasing threshold for 
the mode with average lincwidth. According to Eq. J5J), 
5v = when sin {L 1 /l amp ) turns to zero at L' /l amp = n. 
This "average" lasing threshold agrees with the diffusive 
lasing threshold derived by Letokhov[l6j. Our calcula- 
tion shows that the (conditional) average mode linewidth 
dv can become smaller than the average mode spacing dv 
before the diffusive lasing threshold is reached. Namely, 
with increasing gain, 5v decreases to dv before reaching 
zero. This means the effective Thouless number 6 can be 
reduced to below 1 by coherent amplification for a system 
that is diffusive in the absence of gain. 

Fig. 13 shows the nonlocal part of spatial intensity CF, 
C(Ar, T^ ip / Tamp ) - \C E (Ar, T^ r /T „ m n )\ 2 ■ , in samples of 
g = 4.4 and 9.0. According to Refs. |l7l.ll8|. spatial vari- 
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FIG. 2: Spectral correlation width 8v as a function of amplifi- 
cation time Tamp (triangles) or absorption time T a bs (squares) . 
8u is normalized to the value without gain or absorption 



= OO), 



Tamp / abs 



to Ta mp . Solid symbols correspond 



to g = 4.4, open symbols to g — 9. The solid curves are 
given by the diffusion theory. The inset schematically shows 
the distribution of the resonant mode linewidth, see text for 
discussion. 
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FIG. 3: C - Ci at Ar = and Av = in absorbing 
(squares) and amplifying (triangles) systems. Solid symbols 
correspond to g = 4.4, open symbols to g = 9. Solid and 
dashed curves are obtained from Eq.lJ^Jl without any fitting 
parameters. The inset compares the dependence of C — C\ 
on Ar with (F(Ar) + l)/2 (thick line). Thin solid, dotted, 
and dashed lines represent passive, absorbing, and amplifying 
systems as in Fig0 
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ation and absorption contribution should factorize. Ac- 
counting for terms up to 1/g 2 we obtained the "negative 
absorption" expression for nonlocal intensity correlation 

El 13 EH 

C(Ar, S = L/l amp )-\C E (Ar, S )\ 2 = (l + F(Ar)) x (3) 
1 2s(2 — cos 2s) — sin 2s 4 sin 2 s 
4gs sin 2 s g 2 s 2 

( 2s 2 - s cot s + 1 3 s 2 + s cot s + 1 3s 2 V 
\ 16 sin 2 s 16 sin 4 s 8 sin 6 s / 

where F(Ar) = \C E (Ar)\ 2 . The inset of Fig. plots 
the profile of C(Ar) — Ci(Ar), normalized to its value at 
Ar = 0. For passive, absorbing and amplifying systems, 
the dependence of C — C\ on Ar are almost the same. 
In particular, the value of C — C\ at Ar — > oo is exactly 
1/2 of that at Ar = 0, in agreement with [1 + F(Ar)] 
dependence. This tells us amplification (absorption) in- 
creases (decreases) the nonlocal correlations at every Ar 
uniformly. Therefore, the enhancement (reduction) can 
be characterized by a single number, e.g., the value of 
C — C\ at Ar = as shown in Fig. [3J In two absorbing 
samples of g — 4.4 and 9, the decrease of nonlocal corre- 
lations is in good agreement with the diffusion prediction. 
For amplifying media, only when the fraction of omitted 
lasing realizations is small, Eq.© adequately describes 
the nonlocal correlations of the transmitted intensity. For 
high gain, we see strong deviations: even after removing 
the lasing realizations, nonlocal correlation still exceeds 
the "negative absorption" prediction [Eq. 0]. The de- 
viation becomes more pronounced as the dimensionless 
conductance decreases. The rapid increase of nonlocal 
correlation with gain is caused by enhanced contribution 
from long trajectories that cross upon themselves. 

Finally, we calculated the spectral correlations of 
transmitted intensities. Fig. 0] reveals the changes of 
C\(Av) and C(Av) — C\(Av) due to gain and absorption. 
The frequency separation Av is normalized to the average 
mode linewidth 5v in the absence of gain or absorption. 
C(Av)~ C\{Av) is also normalized to its value at Av = 0. 
Both C\ and C — C\ are narrowed by gain, and broad- 
ened by absorption. While the narrowing of C\(Av) with 
increasing gain is similar to that of a passive system with 
decreasing g, the narrowing of C(A^)—Ci(A^) is just the 
opposite 8] . When a passive system approaches the lo- 
calization transition, the most conducting channels dom- 
inate the transport and correlations. Thus, after normal- 
izing Av by the correlation width 5v, C(Av) — Ci(Av) is 
widened as g decreases in a passive system, but it is nar- 
rowed with the fixed g and increasing gain. This remark- 
able difference illustrates that in a passive system near 
the localization transition, the modes with larger-than- 
average width dominate the transport; while in a diffusive 
system with gain, the modes with narrower-than-average 
width become dominant. 
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FIG. 4: Frequency dependence of nonlocal contribution to 
C(Av) normalized to its value at Av = 0. The inset shows 
the local contribution Ci(Av) = |Ce(A^)| 2 . System param- 
eters and symbol notations are the same as in Fig. The 
open (solid) arrow shows the direction of increasing gain (ab- 
sorption) . 
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